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We present a measurement of the time-dependent CP asymmetry for the neutral B-meson decay 
into the CP — +1 final state A'sTt^tt", with Ks ■k'^iv~ . We use a sample of approximately 227 
million S-meson pairs recorded at the T(45') resonance with the BABAR detector at the PEP-II 
_B-Factory at SLAC. From an unbinned maximum likelihood fit we extract the mixing-induced CP- 
violation parameter S — 0.72±0.71 ±0.08 and the direct CP-violation parameter C = 0.23±0.52 ± 
0.13, where the first uncertainty is statistical and the second systematic. 

PACS numbers: 13.25.Hw, 12.15.Hh, ll.30.Er 



CP violation effects in decays of B mesons dominated 
by 6 ^ sqq transitions (g = m, d, s), are potentially sen- 
sitive to contributions from physics beyond the standard 
model (SM) fl. The B-factory experiments have ex- 
plored time-dependent CP-violating (CPV) asymmetries, 
occuring due to a phase difference between mixing and 
decay amplitudes, in several such decays [1], including 

B° K+K-K° 0, 0, B° /o(980)i4:° Q and 
B" K^Tr° [i]. Within the SM, the magnitude of the 
CPV asymmetry in these decays is expected to be ap- 
proximately equal to the one in & — s- gs decays, such as 
J/^pK'^ A major goal of the B-factory experi- 
ments is to reduce the experimental uncertainties of these 
measurements and to add more decay modes in order to 
improve the sensitivity to beyond-the-SM effects. 

In this paper we present a measurement of the CPV 
asymmetry in the decay B° K°tt°tt°. The K^tt°tt° 
final state is a CP-even eigenstate, regardless of any res- 
onant substructure In the SM this decay is dom- 
inated by the b sqq weak amplitude, with q = u,d, 
and we expect 5 ~ — sin 2/3 and C ~ Here C and 
S are respectively the magnitudes of CP violation in the 
decay and in the interference between decay and mixing, 
and the angle /3 is defined as /3 = arg{—VcdV*ij/VtdVfl), 
where Vij are the elements of the Cabibbo-Kobayashi- 



Maskawa (CKM) quark-mixing matrix pj|. A possible 
contribution from a tree-level b uus amplitude is dou- 
bly Cabibbo suppressed with respect to the leading glu- 
onic penguin diagram. 

The data used in this analysis were collected with the 
BaBAR detector llSll at the PEP-II asymmetric-energy 
e+e" collider A sample of 226.6 ± 2.5 million 

BB pairs was recorded at the T(45) resonance center- 
of-mass energy ^/s — 10.58 GeV. The BABAR detector is 
described in detail elsewhere [T^]. Charged particles are 
detected and their momenta measured by the combina- 
tion of a silicon vertex tracker (SVT), consisting of five 
layers of double-sided detectors, and a 40-layer central 
drift chamber, both operating in a 1.5 T solenoidal mag- 
netic field. Charged-particle identification is provided by 
measurements of energy loss in the tracking devices and 
by an internally reflecting ring-imaging Cherenkov de- 
tector covering the central region. Photons and electrons 
are detected by an electromagnetic calorimeter (EMC) 
composed of 6580 CsI(Tl) crystals. The typical resolu- 
tion for the TT*^ signal in the 77 invariant mass spectrum 



is better than 7 MeV/c^. 

We search for P" K^^tt'^tt^ decays in neutral B me- 
son candidates selected using charged-particle multiplic- 
ity and event topology [l^. We reconstruct tt+tt" 
candidates from pairs of oppositely charged tracks. The 
two-track combinations must form a vertex with a 
probability greater than 0.001 and a 7r"'"7r~ invariant mass 
within 11.2 MeV/c^ of the nominal K° mass [ll]. We 
form tt'^ —^ 77 candidates from pairs of photon candi- 
dates in the EMC, where each photon is isolated from 
any charged track, carries a minimum energy of 30 MeV, 
and has the expected lateral shower shape. B meson can- 
didates are formed from K'^tt^tt'^ combinations and con- 
strained to originate from the e'^e" interaction region 
using a geometric fit. We require that the probability 
of the fit, which has one degree of freedom, is greater 
than 0.001. We extract the K'^ decay length Lj^o and 
the tt'^ — !■ 77 invariant mass from this fit and require 110 
< rrij^ < 160 MeV/c^ and P^o greater than five times 
its uncertainty. The cosine of the angle between the di- 
rection of the decay photons in the center-of-mass system 
of the mother 7r° and the 7r° flight direction in the lab 
frame must be less than 0.92. 

We reconstruct a P° decaying into the CP eigenstate 
ifs7r°7r° (Pcp) and the vertex and flavor of the other 
B meson (Ptag)- The difference At = tcp — ttag of the 
proper decay times is obtained from the measured dis- 
tance between the Bcp and Ptag decay vertices and from 
the boost {l3j = 0.56) of the e+e" system. Ignoring res- 
olution effects, the At distribution is given by: 



V±{At) = 



;-|At|/rgo 

4tro 



-[1=F Aw± 



(1) 



(1 - 2w) (S'sin(ATOdAt) - Ccos(AmdAt))]. 

The upper (lower) sign denotes a decay accompanied 
by a P*^ {B'^) tag, t^o is the mean neutral P lifetime, 
Amrf is the mixing frequency, and the mistag parameters 
w and Aw are the average and difference, respectively, of 
the probabilities that a true P" is incorrectly tagged as a 
P° and vice versa. The tagging algorithm [l7| has seven 
mutually exclusive tagging categories of differing purities 
including one for untagged events that we retain only 
for yield determinations. The effective tagging efficiency, 
defined as the tagging efficiency times (1 — summed 
over all categories, is (30.5 ± 0.6)%, as determined from 
a large sample of P°-decays to fully-reconstructed flavor 
eigenstates (Pflav)- 
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We use the same technique developed for K'^^'^ 
decays of Ref. [§| to reconstruct the K^tt^tt^ vertex 

using the knowledge of the K'^ trajectory and the average 
interaction point (IP) in a geometric fit. The extraction 
of At has been extensively validated in data 9], and on 
large samples of simulated B'^ K'^tt'^tt'^ decays with 
different values of S and C. 

The per-event estimate of the uncertainty on At, 
a{At), reflects the strong dependence of the At resolution 
on the K'^ flight direction and on the number of SVT lay- 
ers traversed by the decay daughters. In about 70% 
of the events both pion tracks are reconstructed from at 
least 4 SVT hits, leading to sufficient resolution for the 
time-dependent measurement. The average At resolu- 
tion, a At, in these events is about 1.0 ps. For events that 
have fewer than 4 SVT hits or for which a{At) > 2.5 ps 
or At > 20 ps, the At information is not used. However, 
since C can also be extracted from flavor tagging infor- 
mation alone, these events still contribute to the mea- 
surement of C. 

We extract the signal yield, S and C from an unbinned 
extended maximum likelihood fit where we parameter- 
ize the distributions of several kinematic and topological 
variables for signal and background events in terms of 
probability density functions (PDFs). 

For each B meson candidate we compute two kine- 
matic variables, the energy difference AE — Eg — 
iy^s and the beam-energy-substituted mass toes = 

^ {^s + po ■ pb)'^/Eq — p\ 13], where the subscripts 

and B refer to the initial T{AS) and the Bcp candidate 
in the lab frame, respectively, and the asterisk denotes 
the e+e^ center-of-mass frame. For signal events, AE 
is expected to peak at zero and mEs at the known B 
meson mass. From a detailed simulation we expect a sig- 
nal resolution of about 3.6 MeV/c^ in mss and 45 MeV 
in AE. Both distributions exhibit a low-side tail due 
to the response of the EMC to photons. We remove a 
small dependence of the signal AE resolution on the lo- 
cation in the X°7r°7r° Dalitz plot by using AE/a{AE) as 
a discriminating variable instead of AE, where a{AE) 
is the calculated uncertainty in AE. We select candi- 
dates with TOES > 5.20 GeV/c^ and -5 < AE/a{AE) < 
2. To suppress other B meson decays we also require 
—0.25 < AE < 0.1 GeV, which does not affect the signal 
AE/a{AE) distribution. 

The background B meson candidates come primarily 
from random combinations of and neutral pions pro- 
duced in continuum events of the type e+e" qq, where 
q — u,d, s, c. Background from BB events may occur ei- 
ther in charmless decays with a K'^ as a decay product, 
or from decays where the Kg is from an intermediate 
charmed particle. The shapes of event variable distri- 
butions are obtained from signal and background Monte 
Carlo (MC) samples and high statistics data control sam- 
ples. The charmless B background forms a broad peak in 



Toes near the S-meson mass; other B background distri- 
butions do not peak in toes- None of the B backgrounds 
peak in the AE / a{AE) distribution. 

We reduce continuum background events, while retain- 
ing 90% of the signal, by requiring | cos^tI < 0.9, where 
9t is the angle between the thrust axis of the Be p can- 
didate's decay products and the thrust axis formed from 
the other particles in the event. We combine 6t, the 
angle between the Bcp momentum and the beam axis, 
9b, and the sum of the momenta pi of the other parti- 
cles in the event weighted by the Legendre polynomials 
Lo{cos{9i)) and L2{cos{9i)) in a neural network (NN). 
The NN has two hidden layers with 4 neurons each, and 
is trained and evaluated [19| on different subsets of sim- 
ulated signal and continuum events and on data taken 
about 40 MeV below the nominal center-of-mass energy. 
To parameterize the NN shape, we divide the NN out- 
put into intervals, chosen such that they are uniformly 
populated by signal events (see, e.g., Ref. [20|). 

We suppress background from other B decays by ex- 



> 



eluding several invariant mass intervals: m{K^Tr 
4.8 GeV/c^ eliminates B° K^tt^ , 1.75 < m{KlTT^) < 
1.99 GeV/c^ reduces B° ^ D°tt° to fewer than 10 ex- 
pected candidates, ■m{n°n^) < 0.6 GeV/c^ removes rjK'^ 
and r]'K°, and 3.2 < TO(7r°7r° ) < 3.5 GeV/c^ removes 
Xco^s and Xc2Kg candidates. 

The signal reconstruction efficiency after all of the 
above requirements is about 15%. Based on MC sim- 
ulations we expect more than one Bcp candidate in 13% 
of the signal events. The selection of the best candidate is 
based only on 7r° information, since the number of mul- 
tiple iCg candidates is negligible (less than 0.1%). We 
select the candidate whose two 7r°s have masses that are 
closest to the expected value. 

For each selected B° K^n'^Tr'^ candidate we exam- 
ine the remaining tracks in the event to determine the 
decay vertex position [l^ and the flavor of iJtag III- 
We parameterize the performance of the tagging algo- 
rithm in a data sample of fully reconstructed i?'^ 
£)(*)-7]-+/p+/u+ decays. For the continuum background, 
the fraction of events in each tagging category is ex- 
tracted from a fit to the data. 

By exploiting regions in data that are dominated by 
background, and by using simulated events, we verify 
that the observables are sufficiently independent that 
we can construct the likelihood from the product of 
one-dimensional PDFs, apart from the signal toes and 
AE /(j(AE) which are correlated. For these observables, 
we use a two-dimensional PDF derived from a smoothed, 
simulated distribution. We obtain the PDF for the At of 
signal events from the convolution of Eq.((T]) with a res- 
olution function TZ(St = At — Attme, CAi)) where A^truo 
is the actual At in the simulated event. The resolution 
function is parameterized as the sum of two Gaussians 
with a width proportional to the reconstructed a At , and 
a third Gaussian with a fixed width of 8 ps fl5| . The first 
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two Gaussian distributions have a non-zero mean, pro- 
portional to (TAt, to account for a bias induced by charm 
decays on the i?tag side. We have verified in simula- 
tions that the parameters of Tl{St, (Ta*) for if^Tr^Tr^ 
events are similar to those obtained from the i?flav sam- 
ple, even though the distributions of cta* differ consider- 
ably. We therefore extract these parameters from a fit 
to the iJflav sample. We also use this resolution function 
for the description of background from other charmless B 
decays. While the resolution functions for B decays into 
open charm final states and continuum have the same 
functional form as used for signal events, the parameters 
for the At PDF of the open-charm background are de- 
termined from MC simulation and they are varied in the 
fit to data for the continuum. 

We subdivide the data into the tagging categories fc, 
events with and without At information (sets I and II), 
and those events located in the inside or outside re- 
gion of the Dalitz plot (in and out). The last subdi- 
vision accounts for the higher contribution and differ- 
ent characteristics of continuum background near the 
Dalitz plot boundary. We define the quantity S — 
mm (m^2 J w^3, 77123)1 where rriij is the invariant mass of 
the Bcp decay daughters i and j combined. This S cor- 
responds to the distance of an event in the Dalitz plot to 
the nearest Dalitz plot boundary in the limit of massless 
daughters. We split the data at 6 = 3.5 GeV^/c^. 

We maximize the logarithm of the extended likelihood 
C = e-(^«+^«) • nLi A- with Ns and Nb = Eb^b 
the total signal and background yields, respectively. The 
likehhood Ck in each tagging category k (with tagging 
fraction e^-) is given as: 

A^I out k 

3 

B 

Nl in k 

n [iV5ef//(l-/L)^A 
j 

ef ^ font) PkAn.j 

B 

Nil out k 

3 

Y.nBe^il-fP)CtQtut,3 

B 

NU in k 

n [iV5ef(l-//)(l-/L)QL- + 
3 

5]nBef(l-/f)(l-/f„,)Q^„,,, 

B 

The probabilities (Q^) and (Q^) for each mea- 



jB 

k,out.j 



jS 



surement j are the products of PDFs for signal (S) and 
background (B) classes: Pk = PDF{mES, AE/a{AE)) ■ 
PDF{NN) ■ PDi^( At, cr( At), tagfc,/c), where for the 
background PDF{mES, AE/a{AE)) = PDF{mEs) ■ 
PDF{AE/<j{AE)). The probabilities Q do not depend 
on At and cr(At) and are used to extract C from the 
yields. The fractions of events with At information for 
signal and background, // and ff, and fractions of 
events in the outside Dalitz plot region, and /J^^ , are 
varied in the fit except for the fractions for B backgrounds 
which are determined from simulation. For about 22% of 
our signal B candidates one or two of the tt'^ decay pho- 
tons associated with Bcp originate from the i?tag- Ac- 
cording to Monte Carlo simulation studies we expect to 
measure the same S and C in these cross-feed events as in 
the correctly reconstructed signal [true) since the contri- 
bution of the 7r° to the At measurement is marginal. To 
account for differences in the PDF distributions for the 
signal probabilities P^ (Q^) we define the signal prob- 
ability to be a linear combination of the correctly re- 
constructed signal and cross-feed events with the relative 
weight determined from simulation. Parameters of signal 
PDFs are the same for the different Dalitz plot regions. 
The PDFs for B backgrounds are identical for the Dalitz 
inside and outside regions. The tagging fractions for the 
signal and the B decay backgrounds are the same, while 
those of the continuum background are different. 

The central values of S and C were hidden until the 
analysis was complete. From a data sample of 33 058 
B° K'^Tr°Tr° candidates, we find Ng = 117 ± 27 
signal decays with S = 0.72 ± 0.71 ± 0.08 and C = 
0.23±0.52 ±0.13 where the first uncertainty is statistical 
and the second systematic. The linear correlation coeffi- 
cient between the two CP parameters is 2%, and the sta- 
tistical significance of the signal yield is 5.8tT. The yield 
of charmless B background is consistent with zero, and 
the fraction of the signal in the outside Dalitz region is 
0.78 ±0.07. Figure [1] shows the distributions of the event 
variables mss, AE/a{AE), and NN output, and the 
ratio of the signal likelihood to signal-plus-background 
likelihood with all variables included. Figure [2] shows the 
At distributions for the 5°- and the S^-tagged subsets, 
and the raw asymmetry [A'^^o — Nbo]/[Nbo + Nbo], where 
the A'^^o {Nbo) is the number of i?'^ (i3°) -tagged events. 
In all plots, data are displayed together with the result 
from the fit after applying a requirement on the ratio 
of signal likelihood to signal-plus-background likelihood 
(computed without the variable plotted) to reduce the 
background. 

We consider the systematic uncertainties listed in Ta- 
ble HI These include the uncertainties in the parameter- 
ization of PDFs for signal and backgrounds which were 
evaluated by varying parameters within one standard de- 
viation or using alternative shape functions. The largest 
contribution to the uncertainty for C is caused by the 
NN shape for continuum inside the Dalitz plot and for S 



7 




5.27 5.29 

mEs(GeV/c2) 



10 12 

AE/o(A E) 



10 





(c)n 

^^^^ 






1 23456789 10 




NN shape 








■ Continuum MC- 




■ B background ; 




Signal MC 




0.2 0.4 0.6 0.8 1 

Ls/(Ls+ Lb) 



1001 



" 0.8 0.85 0.9 0.95 1 
Ls/(Ls+ Lb) 

FIG. 1: Distribution of the event variables (a) niES, (b) 
AE/a{AE), and (c) A^A'^ output in 10 bins after reconstruc- 
tion and a requirement on the ratio of signal likelihood to 
the signal-plus-background likelihood, calculated without the 
plotted variable. The solid line represents the fit result for 
the total event yield and the dotted line for the total back- 
ground. Plot (d) shows the ratio of the signal likelihood to 
signal-plus-background likelihood with all variables included, 
data (dots) with the fit result superimposed. Plot (e) shows 
the same quantity as (d) close to one and with a linear scale. 



from the 2-D parameterization of toes and AE/a{AE). 
We consider uncertainties in the background fractions 
and CP asymmetry in the charmless B background, the 
parameterization of the At resolution function and the 
vertex finding method, knowledge of the event-by-event 
beam spot position, imprecision in the SVT alignment, 
and the possible interference between the suppressed 
b — > ud amplitude with the favored b —^ IJd amplitude 
for tag-side B-decays [IH. We fix Tgo = 1.532 ps and 
Airid = 0.505 ps~^ and vary them by one standard de- 
viation [l^. We correct for the small fit bias which is 
determined using fits to a large number of simulated ex- 
periments, where signal and backgrounds are mixed to- 
gether in the expected proportions. The uncertainty of 
the fit bias is accounted for as a systematic error. 

We perform several consistency checks, including the 
measurement of the B*^ lifetime; we obtain r^o = 1.25 ± 
0.47 ps. We embed different B background samples from 
Monte-Carlo simulation in the data sample and obtain 
consistent yields and CP parameters from the fit. 

In summary, we measure the CP violating asymmetries 




-2 2 
At (ps) 



FIG. 2: Plots (a) and (b) show the At distributions of B°- 
and B°-tagged B° KgTv°Tv° candidates. The solid lines 
refer to the fit for all events; the dashed lines correspond to 
the total background. Plot (c) shows the raw asymmetry (see 
text). A requirement is applied on the event likelihood to 
suppress background. 



in B^ — + K'^TT^TT^ {Kg TT+TT") decays reconstructed 
from a sample of approximately 227 million BB pairs. 
From an unbinned extended maximum likelihood fit we 
obtain S = 0.72±0.71 ±0.08 and C = 0.23±0.52 ±0.13. 
When we fix the values of —5 to the average s'm20 mea- 
sured in b ^ gs modes, sin2/3 = 0.675 ± 0.026 [2^, and 
C to zero, and re-fit the data sample the negative log- 
likelihood changes by 2.2a. The signal yield is consistent 
with our findings in the i^gTr+Tr" decay [l^l assum- 

ing isospin symmetry and that the dominant charmless 
final states are foi980)K°, K*{892)tt°, K*{U30)tt°, and 
non-resonant K'^tt'^tt'^. 



TABLE I: Sources of systematic uncertainties on S and C. 
The total error is obtained by summing the individual errors 
in quadrature. 



Source 


o{S) 


o{C) 


Signal and background PDF parameterization 


0.05 


0.11 


Background fractions 


0.03 


0.02 


CP in charmless B background 


0.03 


0.01 


Vertex finding/resolution function 


0.02 


0.05 


Beam spot position 


0.00 


0.00 


SVT alignment 


0.02 


0.01 


Tag side interference 


0.00 


0.01 


Anid, TB 


0.02 


0.01 


Fit Bias 


0.04 


0.02 


Total systematic error 


0.08 


0.13 
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